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Abstract 



Q^l The minimal supersymmetric standard model with universal boundary 

1^ . conditions and "asymptotic" Yukawa unification is considered. The full 

one-loop effective potential for radiative electroweak symmetry breaking 
as well as the one-loop corrections to the charged Higgs boson, 6-quark 
and r-lepton masses are included. The CP-even Higgs boson masses are 
corrected to two-loops. The relic abundance of the lightest supersymmet- 
d ' ric particle (bino) is calculated by including its coannihilations with the 

next-to-lightest supersymmetric particle (lightest stau) consistently with 
Yukawa unification. The branching ratio of 6 — > 57 is evaluated by in- 
corporating all the applicable next-to-leading order QCD corrections. The 
bino-stau coannihilations reduce the bino relic abundance below the upper 
bound from cold dark matter considerations in a sizable fraction of the 
parameter space allowed by 6 — > 57 for fi > 0. Thus, the /x > case, which 
also predicts an acceptable 6-quark mass, is perfectly compatible with data. 
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It is well-known that the assumption that all three Yukawa couplings of the third 
family of quarks and leptons unify "asymptotically" (i.e., at the grand unified theory 
(GUT) mass scale Mqut ~ 10^^ GeV) naturally restricts the top quark mass to large 
values compatible with the present experimental data. Such a Yukawa unification can 
be obtained by embedding the minimal supersymmetric standard model (MSSM) in a 
super symmetric (SUSY) GUT with a gauge group such as 5'f/(4)c x SU{2)l x S'f/(2)^, 
50(10) or Eq which contain SU{4)c and SU{2)r. Assuming that the electroweak Higgs 
superfields ifi, H2 and the third family right-handed quark superfields f^, If form SU {2)r 
doublets, we obtain the "asymptotic" Yukawa coupling relation ht = hi, and, hence, 
large tan (3 ^ nit/ nib. Moreover, if the third generation quark and lepton SU (2) l doublets 
(singlets) Qs and L3 {If and r'^) form a SU (4)c 4-plet (4-plet) and the electroweak Higgs 
Hi which couples to them is a S'f/(4)c singlet, we obtain hr = hb and the successful 
"asymptotic" mass relation = rrib follows. 

The simplest and most restrictive version of MSSM with gauge coupling unification 
is based on the assumption of radiative electroweak symmetry breaking with universal 
boundary conditions from gravity-mediated soft SUSY breaking. The tantalizing ques- 
tion is then whether this scheme is compatible with exact "asymptotic" unification of 
the three third family Yukawa couplings. A positive answer to this question would be 
very desirable since it would lead to a simple and highly predictive theory. This issue 
has been systematically studied in Ref. 

A significant problem, which may be faced in trying to reconcile Yukawa unification 
and universal boundary conditions, is due to the generation of sizeable SUSY corrections 
to the 6-quark mass @,|]. The sign of these corrections is opposite to the sign of the 
MSSM parameter /i (with the conventions of Ref. 0). As a consequence, for /i < 0, 
the tree-level value of nib, which is predicted from Yukawa unification already near its 
experimental upper bound, receives large positive corrections which drive it well outside 
the allowed range. However, it should be noted that this problem arises in the simplest 
realization of this scheme. In complete models correctly incorporating fermion masses and 
mixing, nib can receive extra corrections which may make it compatible with experiment. 
Also, small GUT threshold corrections to gauge coupling unification can help to reduce 
nib- (For a brief discussion of the possibilities to remedy the nib problem encountered in 
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the ^ < case and some relevant references see Ref. |^.) So, we do not consider this 
6-quark mass problem absolutely fatal for the < case. 

Be that as it may, it is certainly interesting to examine the alternative scenario with 
/i > too. The 6-quark mass receives negative SUSY corrections and can easily be 
compatible with data in this case. This scheme, however, is severely restricted by the 
recent experimental results on the inclusive decay 6 ^ 57 0. It is well-known that 
the SUSY corrections to the inclusive branching ratio BR(6 — 57), in the case of the 
MSSM with universal boundary conditions, arise mainly from chargino loops and have the 
same sign with the parameter fi. Consequently, these corrections interfere constructively 
with the contribution from the standard model (SM) including an extra electroweak Higgs 
doublet. However, this contribution is already bigger than the experimental upper bound 
on BR(6 — i> 57) for not too large values of the CP-odd Higgs boson mass rriA- As a result, 
in the present context with Yukawa unification and hence large tan j3, a lower bound on 
uiA is obtained for /i > 0. On the contrary, for /i < 0, the SUSY corrections to 
BR(6 s'-f) interfere destructively with the SM plus extra Higgs doublet contribution 
yielding, in most cases, no restrictions on the parameters. 

An additional constraint results from the requirement that the relic abundance 
^LSP h"^ of the lightest supersymmetric particle (LSP) in the universe does not exceed 
the upper limit on the cold dark matter (CDM) abundance implied by cosmological 
considerations {^lsp is the present energy density of the LSPs over the critical en- 
ergy density of the universe and h is the present value of the Hubble constant in units 
of 100 km sec~^ Mpc~^). Taking both the currently available cosmological models with 
zero/nonzero cosmological constant, which provide the best fits to all the data, as equally 
plausible alternatives for the composition of the energy density of the universe and ac- 
counting for the observational uncertainties, we obtain the restriction Qlsp h"^ ~ 0.22 
(see Refs. Assuming that all the CDM in the universe is composed of LSPs, we 

further get ^lsp ^ 0.09. 

The LSP is normally the lightest neutralino (x)- In the particular case of Yukawa 
unification, this neutralino turns out to be an almost pure bino. Its relic abundance has 
been estimated, for > 0, and shown to be well above unity, thereby overdosing the 
universe for all m^'s permitted by 6 ^ sj. So, the combination of CDM considerations 
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and the data on BR(6 57) seems to rule out the MSSM with fi > 0, Yukawa unification 
and radiative electroweak breaking with universal boundary conditions. 



It is important to note that, in Refs. the coannihilation |]TT| of the LSP with 

the next-to-lightest supersymmetric particle (NLSP) has been ignored and only the LSP 
annihilation processes have been taken into account. It is well-known, however, that these 
coannihilations can be extremely important, if the mass of the NLSP is relatively close 
to the mass of the LSP, resulting to a considerable reduction of the LSP relic abundance 
||5,12,13]. The question then arises whether, by employing LSP-NLSP coannihilation, one 
can succeed reducing Qlsp below 0.22 for some values of niA allowed by 6 — S7 in 
the ^ > case. This would revitalize a part of the available parameter space for /i > 
saving the simple, elegant and predictive MSSM with universal boundary conditions and 
Yukawa unification even in its simplest realization (with no need of extra corrections to 
nib)- Although, in this parameter range, the sparticles would be quite massive (due to 
the lower bound on niA from b 57) to be of immediate phenomenological interest, we 
would consider this as a very positive development. 

In this paper, we reconsider the constraints from rrib, b ^ s'y and the LSP relic density 
in the context of MSSM with universal boundary conditions and Yukawa unification by 
incorporating the LSP (x) and NLSP coannihilation in the calculation of Qlsp h'^- The 
NLSP turns out to be the lightest stau mass eigenstate f2 and its coannihilation with 
bino has been studied in Ref. (or ||13|) for large (or small) tan/3. Although our analysis 



covers both signs of yU, our main interest here is to see whether, for > 0, there exists a 
range of parameters where all these constraints are simultaneously satisfied. 

We will consider the MSSM with Yukawa unification described in detail in Ref. 
and closely follow the notation as well as the renormalization group (RG) and radiative 
electroweak symmetry breaking analysis of this reference. The only essential improvement 
here is the inclusion of the full one-loop radiative corrections to the effective potential for 
the electroweak symmetry breaking which have been evaluated in Ref. (Appendix 
E). We also incorporate the one-loop corrections to certain particle masses from the 
same reference and the two-loop corrections to the CP-even neutral Higgs boson masses 
(see below). In Ref. 0, we used a constant common SUSY threshold Ms = 1 TeV, 
where the RG-improved tree-level potential was minimized and the parameters n and 
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niA were evaluated. Also, the MSSM RG equations were replaced by the SM ones below 



m 



i2' 



the same scale M5. Here, we use a variable common SUSY threshold Ms = y/mf^ 
where ii^2 are the two stop mass eigenstates. As it turns out, this does not make a very- 
significant numerical diff'erence since the values of /i and niA are found to be pretty stable 
to variations of M5. Despite this fact, it is more appropriate to use a variable SUSY 
threshold here because of the wide variation of the sparticle spectrum encountered. This 
variation appears since, in this work, we consider not only relatively small but also quite 
large values of rriA- Finally, note that our choice of the SUSY threshold Ms minimizes the 
size of the one-loop corrections to fi and ttla and, thus, the errors in the determination 
of these parameters. 

We take universal soft SUSY breaking terms at Mqut, i-e., a common mass for all 
scalar fields mo, a common gaugino mass M1/2 and a common trilinear scalar coupling 
^0, which we put equal to zero (we will discuss later the infiuence of non-zero Aq's). Our 
effective theory below Mqut then depends on the parameters (/iq = h^Mqut)) 

mo, Ml/2, /^o, «G, Mgut, ho, tan/? , 

where ac = Qg/^'^ idc being the GUT gauge coupling constant) and Hq is the common 
top, bottom and tau Yukawa coupling constant at Mqut- The values of ac and Mqut 
are obtained as described in Ref. 0. 

It was pointed out that, for every m^, the requirement of successful radiative 
electroweak symmetry breaking with Yukawa unification implies a relation among the 
GUT parameters tuq and Mi/2- This relation combined with any given ratio of the masses 
of the two lightest SUSY particles (lightest neutralino and stau) leads to a complete 
determination of the values of mo and M1/2 and hence the whole SUSY spectrum for every 
value of niA (see also Ref. [|1^]). We will thus use m^ and = (m^j — m^)/m^ as our 



basic independent parameters. For reasons to become obvious later, we concentrate here 
on the limiting case with m^j = m^ (Afj = 0) where the LSP and NLSP coannihilation 
is [0 most efficient. The values of mo and M1/2 can then be found as functions of uia 
and are depicted in Fig.|I|, for /i > 0, together with the LSP mass {rn^) and the SUSY 
threshold mass parameter Ms- Note that these mass parameters are affected very little 
by changing the sign of /i. 



4 



We find that the effect of including the full one-loop effective scalar potential from 
Ref. is not very significant. In particular, the correction to the tree-level value of 



is found to range from about —7.5% to about 1.5% for between 100 and 700 GeV 
with only a small dependence on the sign of /i. Also, the one-loop radiative correction to 
the tree-level value of fi turns out to be of the same sign as fi and less than about 3.5% 
in the same range of rriA. 

The values of the unified Yukawa coupling constant Hq at Mqut and of tan (3 at Ms 
are estimated by using the running top quark mass at mj, mt{'mt) = 166 GeV, and 
the running tau lepton mass at mz, mr{mz) = 1.746 GeV. We also incorporate the 
SUSY threshold correction to mr{Ms) from the approximate formula of Ref. |]I4|. This 
correction arises mainly from chargino/tau sneutrino loops, is almost m^i-independent 
and has the same sign as /i. It is about 8%, for /i > 0, leading to a value of tan /? = 
55.4 — 54.5 for ttia = 100 — 700 GeV, while, for yU < 0, we find a correction of about —7% 
and tan/3 = 47.8 — 46.9 in the same range of rriA- 

The tree-level values which we find for mb{mz) are quite close to its experimental 



upper bound |T6| : 



mb(m^) = 2.67 ±0.50 GeV. 

The SUSY correction to the bottom quark mass is known to be very large for 
models with Yukawa unification. This correction originates mainly from squark/gluino 
and squark/chargino loops and has sign opposite to the one of (in our convention). 
Thus, for /i < 0, the corrected mf,(m^) will certainly be outside the experimentally 
allowed range. For /i > 0, however, this large negative correction may easily make 
fTT-bif^z) compatible with data. Indeed, using the approximate formula of Ref. [|I4| 



we 



find that, for /i < 0, the correction is about 27.2% - 23.8% for = 100 - 700 GeV, 
which added to a tree-level value of mi,{mz) ~ 3.41 GeV leads to an unacceptably large 
miyimz). For > 0, the tree- level value of mb{mz) is equal to about 3.13 GeV and the 
SUSY correction, to be subtracted from it, about 26.6% — 24.7% for = 100 — 700 GeV. 
The resulting bottom quark mass is then perfectly acceptable in this case. Note that 
the variation of the tree-level value of mb{mz) with the sign of /i is due to the SUSY 
corrections to rrir considered in calculating tan (3. 
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We incorporate in our calculation the two-loop corrections to the CP-even neutral 
Higgs boson mass matrix by employing the program FeynHiggsFast ||T^. The tree- level 
mass of the lightest neutral CP-even Higgs boson is very close to mz for all m^'s consid- 
ered here. Two-loop corrections, however, increase its value, rrih, considerably. The lower 
experimental bound on m^, say equal to 105 GeV, corresponds to rriA ~ 104 GeV. The 
value of rrih increases rapidly with reaching about 115 GeV at = 120 GeV. After 
this, the growth of rrih slows down drastically and this mass soon enters into a plateau 
with rrih ~ 122 GeV. The difference between the two-loop and the tree-level value of the 
heavier CP-even neutral Higgs boson mass mn is insignificant for values of itla bigger 
than about 120 GeV. For smaller m^'s, this difference increases reaching about 10% 
of the tree-level value of rriH at itla ~ 100 GeV. Finally, we also include the one-loop 
corrections to the charged Higgs boson mass using the formalism described in the appen- 



dices of Ref. but without the neutralino and chargino contributions since, as it turns 
out, these contributions are relatively unstable as we change the scale Ms- Presumably, 
higher order corrections will alleviate this instability. It is, certainly, reassuring for our 
procedure that similar results are obtained by using the approximate formula of Ref. [|18|. 
We find that one-loop corrections increase the tree-level value of the charged Higgs boson 
mass mH+- This increase ranges from 7% to 5% of the tree-level value as mA changes 
from 100 to 700 GeV. 

To study the constraints imposed by 6 — 57 on the parameter space of our model, 
we follow the analysis of Ref. [l^. We consider the SM contribution to the inclusive 
branching ratio BR(6 57) from a loop with a W^-boson and top quark (t), the con- 
tribution from loops with charged Higgs bosons (with mass corrected to one-loop) and t 
and the dominant SUSY contribution arising from loops with charginos and stop quarks. 
The SM contribution, which is factorized out in the formalism of Ref. |jT9|, includes the 
next-to-leading order (NLO) QCD |2| and the leading order (LO) QED [||,|l|] 



correc- 



tions. The NLO QCD corrections to the charged Higgs boson contribution are taken 



from the first paper in Ref. ||22|. The SUSY contribution is evaluated by including only 
the LO QCQ corrections p|j23[1 using the formulae in Ref. NLO QCD corrections to 
the SUSY contribution have also been discussed in Ref. [^, but only under certain very 
restrictive conditions which never hold in our case since the lightest stop quark mass is 
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comparable to the masses of the other squarks and the gluinos. Also, the charginos are 
pretty heavy. We, thus, do not include these corrections in our calculation. (Moreover, as 
pointed out in Ref. the available NLO QCD corrections to the SUSY contribution 
are not applicable to models with large values of tan/?, which is our case here.) The 
results, evaluated with central values of the input parameters and the renormalization 
and matching scales, are depicted in Fig.|| for both signs of fi, Aq = and rrif^ = (see 
below). The charged Higgs contribution, added to the one of the SM, raises the predicted 
value of BR(6 57) above the experimental upper bound 4.5 x 10~^ |^ for not too large 
values of m^- The SUSY contribution, which becomes less important as itla increases, 
interferes constructively or destructively with the other two contributions for ^ positive 
or negative respectively. 

We see from Fig.|| that the SM plus charged Higgs contribution decreases as uia 
(or mH+) increases and enters into the experimentally allowed range at rriA ~ 295 GeV 
corresponding to mH+ ~ 318 GeV and Ms ~ 1851 GeV. For /i < 0, inclusion of 
the SUSY contribution makes the BR(6 — > 57) compatible with data for all values of 
tua explored here and no useful restrictions on the parameter space are obtained for 
777.^2 — ^^"^ ^0 = 0. (We may, though, obtain restrictions for m^j's higher than 
and/or Aq 7^ 0.) For > 0, however, the SUSY contribution increases the discrepancy 
between the predicted value of this branching ratio and the data. The upper experimental 
limit on BR(6 — > 57) (^ 4.5 x 10~^) is reached at ^ 385 GeV corresponding to 

^ 694 GeV, mo ^ 781 GeV, M1/2 ^ 1512 GeV and Ms ^ 2418 GeV. We conclude 
that, for /i > 0, = and rrir^ = rriy., uia should be greater than about 385 GeV for 
satisfying the constraints from the b —>■ s'-f process. It is important to observe that, for 
values of rrir^ higher than rriy., charginos become even heavier and the lower bound on 
tua decreases slightly but it can never become smaller than the bound (295 GeV) from 
the SM plus charged Higgs contribution. 

The lower bound on rriA found by using central values of the input parameters can 
be considerably reduced if the theoretical uncertainties entering into the calculation of 
BR(6 — > S'-f) are taken into account. These uncertainties originating from the experi- 
mental errors in the input parameters and the ambiguities in the renormalization and 
matching scales are known to be quite significant. The SM contribution alone, which is 
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factorized out, generates an uncertainty of about ±10%. The charged Higgs and SUSY 
contributions can only increase this uncertainty. In particular, the SUSY prediction for 
/i > should not be a line in Fig.^ but rather a band with the appropriate error margin. 
Consequently, it should intersect the upper experimental bound line at a segment rather 
than at a point. This segment is found to be from about 300 to 510 GeV if only the am- 
biguities from the SM contribution are taken into account. We see that the lower bound 
on rriA is reduced from about 385 to about 300 GeV at least. Inclusion of the errors 
from the charged Higgs and chargino contributions can reduce this bound even further. 
However, we believe that they cannot be reliably calculated at the moment since the 
NLO QCD corrections to the SUSY contribution are not known in our case. In any case, 
the lower bound of 300 GeV on rriA is more than adequate for our purpose here, which 
is to revitalize the MSSM with Yukawa unification, universal boundary conditions and 
fi > 0. Note that, for values of rrif^ higher than rriy. this lower bound on rriA decreases 
slightly but it can never become smaller than the corresponding bound (~ 200 GeV) 
derived from the SM plus charged Higgs contribution. 

The relic abundance of the LSP (x), which is a nearly pure bino in our model, can 
be calculated by employing the analysis of Ref. which is appropriate for Yukawa 
unification and, thus, large tan f3. The inclusion of coannihilation effects of the LSP with 
the NLSP, which is the lightest stau mass eigenstate (f2), is of crucial importance. These 
effects can reduce considerably the LSP relic abundance ^lsp h"^ so as to have a chance 
to satisfy the upper bound 0.22, derived from CDM considerations, for values of rriA 
consistent with the constraints from b S'j. In order to achieve maximal coannihilation 
and thus obtain the strongest possible reduction of Qlsp h"^, we consider the limiting 
case = (see Ref. [^]). The relevant coannihilation processes and Feynman graphs 
together with the corresponding analytical expressions can be found in Ref. [^. Our 
computation here is, however, more accurate since it includes the one-loop corrections to 
the parameter fi, the r-lepton mass, ttla, mH+ and the two- loop corrections to the CP- 
even neutral Higgs boson masses. The results are shown in Fig.|^ for /i > and Aq = 0. 
For comparison, we also include the value of Qlsp obtained by ignoring coannihilation 
effects. Note that these results remain essentially unaltered by changing the sign of fi. 

From Fig.^, one readily finds that the CDM constraint on the LSP relic density 
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0.09 < ^LSP h"^ ~ 0.22 is satisfied for m^'s between about 275 and 400 GeV for /i 
positive, Aq = and nif^ = m^. Of course, the upper bound on uia is more general than 
the lower one since it hold even if the CDM does not solely consist of LSPs. Note that 
the upper bound on (~ 400 GeV) in the /i > case corresponds to m^^ ~ 724 GeV, 
mo ^ 815 GeV, M1/2 ^ 1575 GeV and Ms ^ 2513 GeV. 

It is obvious, from Fig.^, that the reduction of ^lsp h"^ caused by the coannihilation 
effects is dramatic and can bring the LSP relic abundance below 0.22 for m^'s in the 
allowed range from 6 ^ 57 considerations in the fi > case with Aq = and rrif^ = m^. 
Indeed, at rriA ~ 300 GeV, which is its lower bound if a 10% theoretical error is allowed 
in BR(6 57), ^Ilsp h'^ ~ 0.112 (or 3.92 with no coannihilation) and increases with 
niA reaching 0.22 (or 7.4 with no coannihilation) at ~ 400 GeV, where (90% of) the 
central value of BR(6 -> 57) is about (4 x 10""^) 4.44 x 10"^. We see that, for A^^ = 0, 
^0 = and /i > 0, there exists a range of uia (at least between 300 and 400 GeV) 
where both the constraints from b s'y and CDM can be satisfied. Note that, even with 
the central value of BR(6 —>■ 57), this range does not disappear. It only shrinks to the 
interval between 385 and 400 GeV. The value of ^Ilsp with (without) coannihilation 
at TUA ~ 385 GeV is about 0.205 (6.85). 

So far we concentrated in the limiting case = where the coannihilation effects are 
more efficient and we took, for simplicity, Aq = 0. We will now briefiy discuss the effect 
of allowing general values of these quantities. Obviously, for any given tjia, the sparticles 
become heavier as we increase Af^ . The main effect of this is that coannihilation quickly 
faints away and the upper bound on itla from CDM considerations rapidly decreases. 
As a consequence, there exists an upper bound on the parameter Af^ beyond which the 
allowed range of rriA disappears. Positive values of Aq lead to heavier sparticle masses 
and, thus, to an increase of Qlsp h"^ and a slight decrease of BR(6 37). The allowed 
range of ttia again disappears above a positive value of Aq. Negative Aq's bigger than 
about — O.5M1/2 (generally) produce an insignificant decrease in the sparticle masses and 
^LSP h"^- Lower negative values of Aq, however, lead again to an increase of the sparticle 
masses and ^Ilsp h"^, which means that a negative lower bound on Aq must also exist. 

We will not undertake here the difficult task of constructing the region in the mA, 
Af2, Aq space which is consistent with the constraints from b s'j and CDM consid- 
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erations. This would require not only a detailed study of the theoretical uncertainties 
in the calculation of BR(6 — > 57) , but also inclusion of the uncertainties associated with 
the particular implementation of the radiative electroweak symmetry breaking, the RG 
analysis and the radiative corrections to various particle masses. These uncertainties 
propagated to the sparticle spectrum, Qlsp h"^ and BR(6 — >■ 57) can only widen the 
allowed region in the parameter space. Our main conclusion, which is the viability of 
MSSM with Yukawa unification, universal boundary conditions and > 0, can only be 
further strengthened by including these errors. 

In summary, we have considered the MSSM based on radiative electroweak symmetry 
breaking with universal boundary conditions and assumed unification of all three third 
family Yukawa couphngs at the GUT scale. We employed the full one-loop effective 
potential for electroweak symmetry breaking as well as the one-loop corrections to the 
charged Higgs boson, 6-quark and r-lepton masses. Also, two-loop corrections to the CP- 
even Higgs boson masses were taken. We imposed the constraints from 6 — >^ 57 and CDM 
in the universe by carefully including in the LSP relic abundance calculation the LSP 
(bino) and NLSP (lightest stau) coannihilation effects in the case of Yukawa unification 
(and, thus, large tan/9). The calculation of the branching ratio of 6 ^ 57 incorporates 
all the applicable NLO QCD and LO QED corrections and some of its theoretical errors 
were taken into account. We found that bino-stau coannihilation drastically reduces 
the LSP relic density and succeeds to bring it below the CDM upper bound for m^'s 
which are allowed by 6 ^ 57 in the > case. This, combined with the fact that, 
for // > 0, the bottom quark mass after SUSY corrections is experimentally acceptable, 
shows that the simple, elegant and restrictive version of MSSM with Yukawa unification 
and universal boundary conditions can be perfectly viable. It is important to note that, 
without bino-stau coannihilation, this model was excluded. 
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FIG. 1. The values of m^, mo, M1/2 and Mg as functions of niA for > 0, ^0 = and 



= rny. These values are affected very little by changing the sign of 
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FIG. 2. The central value of the SUSY inclusive BR(6 57) as function of for both 
signs of ^, Aq = and ruf^ = fn^. The contributions from the SM and the SM plus charged 
Higgs boson (SM+Higgs) as well as the experimental bounds on BR(6 — 57), 2 x 10~^ and 
4.5 X 10~^, are also indicated. 
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FIG. 3. The LSP relic abundance CIlsp as function of iua in the hmiting case mfj = ruj^ 
and for ij, > 0, Aq = 0. The solid line includes coannihilation of f2 and x, while the dashed 
line is obtained by only considering the LSP annihilation processes. These results are affected 
very little by changing the sign of /x. The limiting lines at 0,lsp h'^ = 0-09 and 0.22 are also 
included. 



